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ABSTRACT 

We  present  a  study  of  the  interaction  of  small  amplitude,  unsteady,  freestream  distur¬ 
bances  with  a  shock  wave  induced  by  a  wedge  in  supersonic  flow.  These  disturbances  may 
be  acoustic  waves,  vorticity  waves,  or  entropy  waves  (or  indeed  a  combination  of  all  three). 
Their  interactions  then  generate  behind  the  shock  disturbances  of  all  three  classes,  an  aspect 
that  is  investigated  in  some  detail,  our  motivation  here  being  to  investigate  possible  mech¬ 
anisms  for  boundary- layer  receptivity,  caused  through  the  amplification  and  modification  of 
freestream  turbulence  through  the  shock-body  coupling.  Also,  the  possibility  of  enhanced 
mixing  owing  to  additional  vorticity  produced  by  the  shock-body  coupling  is  investigated. 
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1.  Introduction. 

Boundary-layer  receptivity  ( i.e  the  growth  of  instabilities  within  the  boundary  layer, 
triggered  by  some  form  of  external  disturbances)  has  been  observed  in  supersonic  wind  tun¬ 
nel  experiments  in  the  presence  of  sound  waves  radiated  by  turbulent  tunnel- wall  boundary 
layers  or  by  freestream  turbulence.  In  such  situations,  the  acoustic  wave  mainstream  dis¬ 
turbances  are  first  processed  by  the  shock  wave  induced  by  the  body  in  supersonic  stream, 
and  then  interact  with  the  boundary  layer  to  cause  instability  (if  at  all).  In  this  paper, 
we  study  the  former  phase  of  how  a  shock  in  the  presence  of  the  body  processes  main¬ 
stream  disturbances.  The  latter  phase  of  these  disturbances  internalizing  in  the  boundary 
layer  will  be  subject  of  future  study.  The  shock-body  combination  is  shown  to  produce 
a  wide  spectrum  of  wavelengths  from  monochromatic  waves,  i.e  a  freestream  disturbance 
with  a  single  wavelength;  in  particular,  we  note  the  production  of  shorter  wavelength 
disturbances  than  would  be  expected  from  the  free  shock-disturbance  interaction.  Then 
these  shorter  wavelength  perturbations  resulting  from  the  shock- body  combination  would 
go  on  to  interact  with  the  boundary  layer  on  the  surface  of  the  body,  leading  to  possible 
instability.  In  the  case  of  turbulence-shock-body  interaction,  it  is  noted  that  the  turbu¬ 
lence  itself  is  composed  of  a  continuous  spectrum  of  wavelengths  and  the  interaction  would 
produce  even  more  short-wavelength  disturbances  than  the  single  wavelength  disturbance 
considered  here.  Such  interactions  could  possibly  enhance  mixing  of  a  multi-component 
flow  owing  to  the  amplification  and  generation  of  vorticity  behind  the  shock.  Again,  it 
is  found  that  the  shock-body  combination  efficiently  processes  the  freestream  disturbance 
and  generates  a  vorticity  profile  behind  the  shock  which  are  larger  than  those  arising  from 
the  interaction  of  the  free  shock  with  mainstream  disturbances. 

To  investigate  these  phenomena,  we  consider  the  model  problem  of  a  wedge  in  a 
uniform  supersonic  stream.  In  the  presence  of  the  wedge,  the  supersonic  flow  abruptly 
changes  direction  through  an  oblique  shock  wave,  details  of  which  are  determined  by  the 
Rankine-Hugoniot  conditions.  This  paper  is  concerned  with  the  response  of  the  overall 
flow  field  to  general,  small  amplitude  disturbances  in  the  freestream  that  are  convected 
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through  the  shock  wave.  In  order  to  quantify  this  response,  we  consider  isolated  single 
wavelength  disturbances  of  acoustic  type,  vorticity  type  or  entropy  type.  In  particular, 
we  are  concerned  with  the  production  of  a  continuous  spectrum  of  wavelengths  produced 
in  the  flow  behind  the  shock  in  the  presence  of  the  wedge.  Moore  (1954),  Ribner  (1953) 
and  McKenzie  and  Westphal  (1968)  have  investigated  the  processing  of  single  wavelength 
disturbances  by  freely  propagating  shocks,  and  note  that  through  the  interaction  of  the 
shock  and  the  freestream  disturbance,  single  wavelength  entropy,  acoustic  and  vorticity 
disturbances  are  all  generally  produced  behind  the  shock  in  response  to  an  isolated  acoustic 
disturbance,  or  to  an  isolated  entropy  disturbance  or  to  an  isolated  vorticity  disturbance. 
With  the  presence  of  the  wedge,  the  flow  behind  the  shock  is  no  longer  comprised  of  single 
wavelength  acoustic,  entropy  and  vorticity  disturbances.  In  order  to  satisfy  the  condition 
of  zero  normal  velocity  at  the  wedge  surface  as  well  as  the  condition  that  the  shock  remain 
attached  to  the  apex  of  the  wedge,  a  continuous  spectrum  of  wavelengths  for  the  acoustic, 
entropy  and  vorticity  disturbances  exists  behind  the  shock.  It  is  this  more  complicated 
disturbance  pattern  that  goes  on  to  interact  with  the  boundary  layer  developing  along  the 
surface  of  the  wedge,  possibly  initiating  new  receptivity  mechanisms. 

Limited  work  appears  to  have  been  done  on  the  theoretical/computational  approach 
to  the  interaction  between  shock  waves  and  boundary-layer  instabilities,  although  most 
of  this  work  has  focused  on  the  wedge  problem.  A  fairly  restrictive  (and  somewhat  ad 
hoc)  model  was  presented  by  Petrov  (1984),  using  inviscid  linear  stability  equations,  and 
somewhat  heuristic  conditions  on  the  shock  (which  was  also  assumed  to  lie  at  the  edge  of 
the  boundary  layer).  Cowley  and  Hall  (1988)  presented  an  asymptotic  model,  applicable 
to  three-dimensional  viscous  modes  of  instability,  with  appropriately  simplified  conditions 
applied  on  the  shock  surface  (derived  from  the  Rankine-Hugoniot  conditions),  which  was 
taken  to  lie  just  outside  of  the  boundary  layer.  Chang,  Malik  and  Hussaini  (1990)  con¬ 
sidered  the  viscous  small  disturbance  equations  with  parallel  mean  flow  assumption,  using 
the  full  inviscid  Rankine-Hugoniot  conditions  on  the  shock. 

Here  we  take  the  equations  governing  the  gas  flow  both  ahead  of  and  behind  the  shock 
to  be  the  Euler  equations  (neglecting  the  boundary  layer  on  the  wedge  in  the  first  order  of 
approximation)  along  with  the  ideal  gas  law.  The  no  mass  flux  condition  is  imposed  on  the 
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wedge  surface,  so  that  the  normal  component  of  the  flow  must  have  the  same  velocity  as 
the  velocity  of  the  wedge  normal  to  its  surface.  Furthermore,  it  is  required  that  the  shock 
remains  attached  to  the  apex  of  the  wedge.  The  analysis  is  accomplished  by  linearizing 
the  Euler  equations  about  the  base  state  and  applying  the  Rankine-Hugoniot  conditions 
at  the  mean  position  of  the  shock  wave. 

Carrier  (1949a)  and  Van  Dyke  (1953)  first  investigated  the  response  of  the  inviscid  flow 
field  to  time-periodic  oscillations  of  the  wedge.  Their  interest  was  in  the  surface  pressure 
distribution  and  the  resultant  forces  and  moments,  and  their  relevance  to  oscillating  airfoils. 
The  related  problem  of  freestream  turbulence  amplification  caused  through  interactions 
with  shock  waves  has  been  considered  by  Anyiwo  and  Bushnell  (1982),  Hussaini,  Collier 
and  Bushnell(1986),  Meadows,  Kumar, and  Hussaini  (1991),  and  Kumar,  Bushnell  and 
Hussaini  (1989)  who  confirmed  the  important  result  found  in  Moore  (1951),  Ribner  (1954) 
and  McKenzie  and  Westphal  (1968)  that  a  pure  acoustic,  vorticity,  or  entropy  wave,  upon 
interacting  with  a  plane  shock  wave,  generally  generates  all  three  classes  of  disturbance 
downstream.  The  interaction  of  a  shear  wave  with  a  detonation  wave  induced  by  a  wedge 
in  a  supersonic  flow,  as  well  as  the  response  of  the  oblique  detonation  to  oscillations  of 
the  wedge  was  considered  by  Lasseigne  and  Hussaini  (1992).  In  studying  the  interaction 
of  a  shear  wave  with  the  detonation,  they  assumed  that  a  weak  steady  sinusoidal  vorticity 
wave  is  obliquely  convected  through  an  overdriven  detonation  attached  to  a  wedge.  Their 
concern  was  with  the  response,  measured  by  the  deviation  of  the  detonation  position  from 
its  unperturbed  state  and  by  the  vorticity  and  pressure  at  the  detonation,  as  the  degree 
of  overdrive  is  increased.  They  compared  the  flow  in  the  presence  of  the  wedge  with  the 
unobstructed  flow  field  considered  in  Jackson,  Kapila,  &  Hussaini  (1990).  In  particular, 
for  this  special  type  of  disturbance  (i.e.  steady),  they  found  an  infinite  discrete  spectrum 
of  disturbances  to  exist  behind  the  shock. 

The  present  study  also  raises  other  important,  broader  issues,  notably  the  stability  of 
the  shock  on  the  wedge.  It  is  generally  accepted  that  the  so-called  “strong  shock”  solution 
(Liepmann  and  Roshko  1957)  is  unstable  (see  Levinson  1945,  Carrier  1949b,  Henderson 
and  Atkinson  1976,  Rusanov  and  Sharakshanae  1980,  and  Salas  and  Morgan  1982),  whilst 
the  proof  for  the  stability  of  the  “weak  shock”  solution  has  largely  been  either  numerical 
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(Rusanov  and  Sliarakshanae  1980,  and  Salas  and  Morgan  1982)  or  subject  to  some  restric¬ 
tions.  The  work  of  Henderson  and  Atkinson  (1976)  considered  just  finite  length  wedges  to 
"avoid  unbounded  velocity  downstream”  (sic)  whilst  Carrier  (1949b)  did  "not  worry  about 
convergence  in  the  large”  (sic)  when  considering  series  solutions.  One  of  the  aims  of  this 
paper  is  to  place  the  evidence  for  the  stability  of  flows  in  which  the  flow  downstream  of 
the  shock  is  supersonic  on  a  much  firmer  footing.  Here  we  shall  remove  Henderson  and 
Atkinsons’s  (1970)  restriction  and  pay  particular  attention  to  the  convergence  of  series 
solutions. 

The  layout  of  the  paper  is  as  follows.  In  Section  2  we  consider  the  three  distinct  classes 
of  disturbances  upstream  of  the  shock,  and  then  in  Section  3,  the  analytic  solution  for  the 
downstream  development  of  these  disturbances  in  the  region  bounded  by  the  shock  and 
the  wedge  surface  is  obtained.  In  Section  4,  we  examine  various  aspects  of  the  analytical 
solution.  These  raise  important  questions  regarding  the  stability  of  the  shock  itself.  We 
show  that  the  weak  shock  solution  is  stable  to  these  imposed  disturbances,  provided  that 
the  downstream  flow  is  supersonic  in  nature,  in  line  with  previous  (mostly  numerical  and 
experimental)  studies.  We  also  examine  the  behaviour  of  the  far- downstream  flow  which 
suggests  an  important  physical  decomposition  of  the  solution.  The  physical  aspects  of  the 
decomposition  is  examined  in  Section  5  and  the  responses  to  the  various  types  of  imposed 
disturbances  are  quantified.  Our  conclusions  are  given  in  Section  6. 


2.  Formulation 

Throughout  we  will  denote  quantities  upstream  of  the  shock  by  subscript  1  and  down¬ 
stream  by  subscript  2.  We  take  the  wedge  surface  to  make  an  angle  8  with  respect  to  the 
on-coming  flow,  with  (x*,y*)  the  coordinates  parallel  and  perpendicular  to  the  upstream 
base  flow  (respectively)  which  has  magnitude  {/,*,  Mach  Number  M\,  density  p *,  and  tem¬ 
perature  T* .  We  assume  that  the  ratio  of  specific  heats  7  is  constant.  The  velocity  vector 
is  written  as  U*{u,v)  with  respect  to  (,r*,y*)  coordinates.  The  density  field  is  written  as 
p*p,  pressure  as  p*R*T*p ,  where  E*  denotes  the  gas  constant  and  the  temperature  field 
as  T*T. 
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The  Rankine-Hugoniot  relations  provides  a  link  between  conditions  upstream  and 
conditions  downstream  of  the  shock,  and  in  particular  leads  to  the  following  classical 
result  in  gas  dynamics  (e.g  Liepmann  and  Roshko  1957), 

tan(/?  -  9)  _  (h  _  (h  _  (7  -  l)^2  sin2  ft +  2 
tan  ft  uj  p2  (7  4-  l)Mf  sin2  ft 

where  ft  is  the  angle  between  the  shock  and  the  wedge  centerline,  and  u\  and  U2  are  the 
non-dimensional  velocity  components  perpendicular  to  the  shock. 

It  is  well  known  (Liepmann  and  Roshko  1957)  that  for  a  given  upstream  Mach  number 
M\  and  9  <  9max(M\)  equation  (2.1)  admits  two  possible  solutions  for  ft:  the  larger  value 
corresponding  to  the  so  called  “strong  shock”  solution,  which  is  characterized  by  subsonic 
flow  behind  the  shock  in  all  cases  and  the  smaller  value  corresponding  to  the  so-called 
weak  solution  which  is  characterized  by  the  flow  being  generally  supersonic  behind  the 
shock,  except  for  a  small  region  of  subsonic  flow  for  9  close  to  9  =  9rnax .  For  9  >  9rnaz , 
no  attached-shock  solutions  to  the  problem  exist.  Attached  shock  solutions  all  have  the 
property  of  uniform  downstream  flow  directed  parallel  to  the  wedge  surface. 

Small  amplitude  disturbances  upstream  of  the  shock  may  be  classified  into  three  dis¬ 
tinct  classes  (see  McKenzie  and  Westphal  1968,  for  example).  Taking  e  to  be  the  measure 
of  the  amplitude,  and  hence  a  small  parameter,  we  have: 

(t)  Acoustic  waves:  These  are  characterized  by  having  a  pressure  perturbation  and  the 
resultant  perturbations  in  velocity,  density  and  temperature;  however,  the  waves  carry 
no  change  in  entropy  and  have  no  vorticity.  In  this  case  the  disturbance  field  upstream 
of  the  shock  is  given  by 


p  =  1  +  eE  -f  0(e2), 

(2.2) 

7A/?(«,+u,)  (  >’ 

(2.3) 

a2eE  2 

TM12(a1+  u)+°{€  )? 

(2.4) 

P  —  Id - -  +  0(e2), 

(2.5) 
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(2.6) 


T=l-  — — —+0(e2), 

1 

where  E  is  the  normal  mode  exponential 

E  —  exp{miori  +  i&2y\  +  (2.7) 

and 

uj  —  —a\  ±  d  (2.8) 

is  the  frequency  of  the  disturbance  for  the  given  wavenumbers.  The  coordinates  j] 
and  t/i  are  parallel  and  perpendicular  to  the  upstream  flow  respectively,  suitably  non- 
dimensionalized  (this  can  be  accomplished  by  using  one  of  the  wavelengths  of  the 
disturbance  as  a  characteristic  length  scale,  e.g.  setting  c*],  say,  to  unity).  The  choice 
of  frequency  u>  and  also  the  nature  of  the  solution  (2.2)-(2.6)  arises  from  solving  the 
irrotational  flow  problem  upstream  of  the  shock.  We  shall  refer  to  modes  correspond¬ 
ing  to  the  negative  sign  in  (2.8)  as  so-called  “fast  modes”,  and  to  the  positive  sign  in 
(2.8)  as  the  so-called  “slow  modes”. 

(ii)  Vorticity  waves:  These  are  characterized  by  having  no  density,  temperature  or  pressure 
components  to  the  disturbance  to  0(e).  As  such,  the  upstream  field  can  be  written 


u  =  l-t-  eE  +  0(e2), 

v=  -e^L  E  +  0(e2), 
«2 

p,T,p=l  +  0(e2), 


together  with 


u>  =  -«!, 


(2.10) 

(2.11) 


(2.12) 


with  E  having  the  same  definition  as  in  (2.7). 


(Hi)  Entropy  waves:  These  are  characterized  by  having  no  pressure  or  velocity  components 
to  the  disturbance  to  0(e).  As  such,  the  upstream  field  can  be  written 


P  =  l  +  eE  +  0(e2), 


(2.13) 
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T  =  l  -  eE  +  0(e2), 

(2.14) 

u,  p  =  1  -f  0(e2), 

(2.15) 

v  —  0(e2), 

(2.16) 

together  with  u  defined  as  in  (2.12)  and  E  defined  as  in  (2.7). 


In  the  following  section  we  go  on  to  consider  the  interaction  between  these  waves 
and  the  shock  attached  to  the  wedge.  Fortunately,  in  spite  of  their  distinct  features,  the 
analysis  is  fundamentally  the  same  in  all  three  cases,  and  wt  shall  see  how  just  one  of 
these  waves  upstream  of  the  shock  generally  produces  a  combination  of  all  three  modes 
of  disturbance  downstream  of  the  shock;  however,  unlike  the  situation  where  no  wedge  is 
present,  the  downstream  disturbance  is  not  restricted  to  a  single  wavelength. 


3.  The  Solution  Downstream  of  the  Shock 

Downstream  of  the  shock,  the  flow  variables  may  be  written 

P  =  Pi  +  cp  +  0(e2),  (3-1) 

p  =  p2  +  ep+ 0(e2),  (3.2) 

T  =  T2  +eT  +  0(e2).  (3.3) 

The  non-dimensional  coordinates  parallel  and  perpendicular  to  the  wedge  are  taken  as  x2 
and  i)2  respectively  whilst  the  velocity  components  in  the  x2  and  y2  direction  are  written 
as  U2  +  eu  4-  0(e2)  and  eu  +  0(e2)  respecti\  3ly.  The  governing  equations  at  0(e)  may  be 
written 


Pt  +  U2PX2  +  P2UXi  +  P2Vy2  —  0, 

(3-4) 

p2{ut  +  U2UX2 }  +  iM2P*.  ~ 

(3.5) 

P2{vt  +  U2VX?}  +  ^  ^2Py2  —  0? 

(3.6) 

p2{Tt  +  U2TI2}-  ^{pt  +  U2pl2)  =  0. 

(37) 

i 


P  -  pzT  +  Tip. 


(3.8) 


TU“  general  solution  procedure  is  based  on  that  followed  by  Carrier  (1949a)  and 
\an  Dyke  (1953).  It  is  found  convenient  to  first  split  the  solution  into  two  components, 
corresponding  to  acoustic  waves  and  vorticity  waves.  Specifically,  we  write 


U  =  <*>*,  +  £  #s:  (3.9) 

e  0y2  El?  i  (3.10) 

where  4>  =  <?>(.r2,  y2,  f )  represents  the  acoustic  inode  and  E  =  E(x2,y2,t)  represents  the 
vorticity  mode. 

Substituting  these  forms  into  (3.5)-(3.6),  differentiating  and  combining  appropriately 
yields  the  following  equations 


V2{*  +  C2«„  +  -^-}  =  0,  (3.11) 

V2{E(4-{J2JSI2}  =0.  (3.12) 

Solving  (3.11)  and  (3.12)  implies  the  introduction  of  two  arbitrary  harmonic  functions  (see 
the  footnote  in  Hui  1969).  Equations  (3.11)  and  (3.12)  taken  together  represent  a  sixth 
order  system  which  is  derived  from  the  second  order  system  (3.5)-(3.6),  and  therefore 
without  loss  of  generality,  we  set  these  two  arbitrary  harmonic  functions  to  zero  and  note 
that  all  physical  conditions  at  the  shock  and  the  wedge  boundary  are  satisfied  by  the 
resulting  solution.  Using  (3.4),  (3.7)  and  (3.S),  it  is  then  possible  to  show  that  <t>  satisfies 
the  following  equation 


with 


where 


—  ~{4>tt  +  2 U2^x2<  +  U 2  $ 1 2X2 } , 


(3.13) 


V 


IPz 


(<f>t  +  lT2<bt ), 


(3.14) 


(l  Cf 


AIi 


(3.15) 
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is  the  nondimensional  sound  speed  behind  the  shock.  A  further  quantity  is  also  required 
to  describe  the  proui^m  fully,  namely  4''{y2,  t),  the  displacement  of  the  shock  wave. 

It  is  possible  to  write  the  general  solution  for  <f>  and  E  which  is  bounded  at  the  tip  of 
the  wedge  as 


where 


lUii 

<j>  =  e 


{fli/  cosh(^02)  +  b„  sinh(iA°2)} 

(3.16) 

i/=0 

(3.17) 

i>=0 

.  I  j  U  2  A  p  o 

0'  =  e  ^  dvJyjkjyz), 

(3.18) 

i/=0 

i fc  =  ■£, 

(3.19) 

al 

(3.20) 

W 

u  = 

a , 

(3.21) 

£=yJv-02, 

(3.22) 

r2  =  xl-02yl 

(3.23) 

tanh  9  2  = 

x2 

(3.24) 

A  —  cot(/?  —  B). 

(3.25) 

In  deriving  (3.16)-(3.18),  we  have  assumed  that  the  harmonic  time  dependency  upstream 
of  the  shock  is  transmitted  unchanged  through  the  shock,  an  assumption  that  is  likely  to  be 
accurate  provided  that  |e|  <<  1  and  the  wavelength  of  the  disturbance  is  large  compared 
to  the  shock  thickness.  The  boundedness  argument  at  the  wedge  tip  (r  =  0)  precludes  the 
Yv(z)  family  of  solutions. 

If  we  impose  the  impermeability  condition  on  the  wedge  surface,  then  the  b v  coeffi¬ 
cients  must  all  be  set  to  zero;  if  we  demand  that  the  shock  remains  attached  to  the  wedge 


9 


tip  then  t/u  =  0.  and  if  we  assume  that  there  is  no  transient  behaviour  (i.c.  only  tin' 
sustained  response  to  the  imposed  disturbance  is  present)  then  «o  =  c0  —  0  also. 

The  key  useful  results  are 


~)p2  t  u><- 

p  =  —e 


iwt/ 

l’i-a 


E'k 

e=0 


9  4.  U> 

2iav  H - at 

a. 


U2 

i 


cosh(  u02  )Jv(kr), 


(3.26) 


iLljt  - 


iuj  x  2  iu/  y  2  Q  j  OO 

“r7r“e2«"l-?.  Y^kl  -iX 


v  =  0 


as  1  -  1  - 

pTCu  d  2‘'rt/+1  —  9^cc-i  r  Jv(k£,y 2) 


1  -<-7^ 

—  —  C  2  * 

9 


^  -  f  .  U 2  \ 

2_^  k  <  2? — u^+ay-]  -  av+\  >  cosh(n02 )Jv{kr), 

t/  =  0  v  J 


(3.27) 


■/y2  Aa*  OO 


r/=0 


ik02as 


cuJu{k(y 2) 


1  iwt-  00 

+  -e  ^-a' 


(3.28) 


^  kfi  +  a^+i  }  sinh(n02)  J„(fcr), 


i/=0 


where  u_i  =  c_  1  —  d-\  =0.  The  set  of  solutions  above  now  involves  three  sets  of  coeffi¬ 
cients,  namely  {«i/ },  and  {d„}.  Although  these  solutions  satisfy  the  impermeability 
constraint  on  02  =  0,  we  have  yet  to  impose  the  unsteady  Rankine-Hugoniot  conditions 
(McKenzie  and  Westphal.  1968)  on  these  perturbation  terms  at  the  shock,  which  serve  to 
determine  the  aforementioned  set  of  coefficients. 

The  Rankine-Hugoniot  conditions  on  the  shock  may  be  written  in  the  following  form 


u'2  +  u2  —  —  ( 1 - )R'  -j - (u\p\  +  «j ), 


P  2 


P2 


p2 


,  ,  !-  Pi  ,  a2s  p2  ,  ,  1_  .  ,  Pi 

2  p2  2  ')u2p2  2  2~iM2u\ 

u'2  _  a*  P±  +  fl*  P2  =  (1  _  +  »i»i 

2  (7  -  1)«2  Pi  (7  -  1)^2  Pi  u2  ii2 


P 1 


+ 


Pi 


(y-l  )A/f  f<2  (7-I )  M2  u  -2  ' 

>’2  =  *'!+(»  1  -  «2)l/’». 


Here  we  have  written 


R'  =  »’i'/v  +  »/v 


(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 
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where  a  denotes  distance  along  the  (unperturbed)  shock.  We  have  taken  ( u',v ')  to  be  the 
velocity  perturbations  perpendicular  and  parallel  to  the  shock,  respectively,  and  (u,v)  the 
basic  flow  perpendicular  and  parallel  to  the  shock  respectively. 

The  four  conditions  (3.2yj-(3.33)  may  be  combined  to  eliminate  the  downstream  per¬ 
turbation  density  p'2.  It  is  then  possible  to  insert  the  expressions  for  p ,  u,v,ip' ,  all  evaluated 
on  the  undisturbed  shock  surface,  into  these  resulting  equations  to  yield  the  following  sys¬ 
tem  of  equations  at  each  value  of  v  >  0. 


cos(/?  —  6) 


{!4r4 


£ci/+i 


w} 


•  1  .U2 

--[2 i  —  a„  —  a^+i  4-  au-\]  cosh^# 

2 

+  sin (/3  -  9)  |  +  a„+i]  sinh(i/0n)  + 

~{ux  -  u2)[-—du  +  1  -  <f„-i)]}  =  R[l\ 

&8  ^  J 

|  |  |sin(/?  -  6) 

—  \[2i— av  —  civ+i  +  a^-i]  cosh(^o)} 

2  a$  J 

-  cos(/3  -  6)  +  a„+i]  sinh(i/#o)  +  j  j 

+  {  ^2  }  sin^  ~  ^  1  -  ^-0]  +  i{d2asdv | 

+  (  I3  I  ( \{2iau  +  —a„- 1  -  —  a^+i]  cosh^flo) 

[  -^3  J  ^  2  flj  ^3 

p(2) 
ill/ 


r  *<2>  i 

1  i?l3)  j 


Here  we  have  written 


tanh  6o  =  0  tan(/?  —  0), 

=  -(i  -  — ), 

2  p-2 

n 

-4.3 


2u2  ’ 


(3.34) 


(3.35a,  b) 


(3.36) 

(3.37) 

(3.38) 

(3.39) 
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(3.40) 


B,  =  1  + 


(7  -  u«r 


-  ((7-l)«l  p2)' 


B,  = 


(7  ~  1)«2 


(3-42) 


The  determination  of  and  is  routine  once  the  upstream  disturbance  is 

prescribed,  but  for  brevity  these  are  omitted.  These  quantities  arise  from  the  forcing  caused 
by  the  upstream  disturbance  which  is  proportional  to  the  exponential  (2.7)  evaluated  at 
the  shock  position.  This  exponential  is  expanded  in  terms  of  the  Bessel  function  by  using 


l/ow  A 

—  *  T7y~TTff2 


exp[i(«i  cot  /3  +  a2)(cos0  +  Asin0)j/2]  =  e  1 2  “*  eiS>n6kl'y2 

'u*—*92  e  +  (-1)e  w;;  \ 

=  e  2  '  2^ - 7X7 - *Mfc£y2), 

^  1  + 


(3.43) 


where 


f  (01  cot  fj  +  o?2)(cos  fl  +  A  sin  6)  4- 
6  =  sin"1  ^  - — - 3—^ 


(3.44) 


and  Si/}q  is  the  Kronecker  delta  function.  Then  a„+i,  6„+i,  and  cv+\  are  determined  from 
equations  (3.34),  (3.35)  which  are  essentially  recursive  in  form. 


4.  Analysis  of  the  Series  Solution  and  Far-Downstream  Behaviour 

In  Figure  la,  we  show  3f?{pe_,w<}  (solid  line)  and  ^s{pe~lu,t }  (broken  line)  both  evalu¬ 
ated  along  the  wedge  surface  for  the  particular  case  M\  =  5,  0  =  25°,  oq  =  1,  and  a 2  =  0, 
taking  the  positive  sign  in  (2.8)  (the  so  called  “slow  mode").  The  corresponding  shock  po¬ 
sition  3?{  ip'c ~,ujt }  (solid  line)  and  (broken  line)  with  distance  along  the  shock, 

r2  —  .1/2(1  +  A2)1/2  is  shown  in  Figure  lb.  These  results  indicate  a  complex  structure  to 
the  solution:  both  quantities  show  the  expected  oscillatory  behaviour,  although  the  wedge- 
surface  pressure  shows  indications  of  amplitude  decay  as  rz  — >  00.  Shown  in  Figure  2  is  the 
response  to  the  so-called  “fast  mode",  i.t.  the  negative  sign  in  equation  (2.8)  is  used.  The 
wedge-surface  pressure  is  shown  in  Figure  2a  whilst  the  shock  position  is  shown  in  Figure 
2b;  as  before,  real  components  are  indicated  by  solid  lines  and  imaginary  components  by 
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broke:;  lines.  The  expected  oscillatory  behaviour  is  shown,  but  in  this  case  there  is  little 
evidence  of  amplitude  modulation  of  the  wedge-surface  pressure.  These  two  examples  indi¬ 
cate  the  importance  of  ascertaining  the  correct  downstream  behaviour  of  the  solution  given 
by  the  infinite  sums  (3.26-3.28).  Unfortunately,  it  was  found  that  as  V2  (or  s?)  became 
large,  the  computations  for  1/  (and  p)  become  extremely  sensitive  to  the  number  of  terms 
of  the  series  taken  (although  we  stress  that  all  results  presented  are  carefully  checked  for 
accuracy,  in  particular  for  insensitivity  to  truncation  of  the  number  of  terms  in  the  Bessel 
series  (3.16)-(3.18)).  Cursory  inspection  of  the  coefficients  au,  b and  c/„  indicated  that 
the  numerical  difficulties  described  result  from  these  coefficients  often  becoming  extremely 
large  as  u  —>  oc.  leading  to  numerical  difficulties  associated  with  the  effects  of  roundoff. 

The  same  general  1 j. end  of  oscillatory,  complex  behaviour  seen  in  Figures  1  and  2  was 
also  seen  in  the  results  of  Lasseigne  and  Hussaini  (1992)  relating  to  detonation  waves  and 
moving  wedges.  However,  there  are  two,  not  unrelated,  issues  that  remain  unanswered  and 
affect  the  physical  interpretation  of  the  results.  These  two  issues  are:  (i)  the  nature  of  the 
coefficients  a „,  c„  and  d„  as  n  — >  oo  and  (ii)  the  precise  nature  of  the  flow  far  downstream 
of  the  shock.  We  look  at  these  aspects  in  some  detail  and  discover  a  number  of  interesting 
subtleties. 

Superficial  inspection  of  the  coefficients  indicates  general  growth  as  u  — ►  oo.  Indeed, 
this  was  observed  even  when  the  inhomogeneous  forcing  terms  R\}\  and  were 
“turned  off’  after  u  =  1,  in  order  to  mimic  a  homogeneous  system  (although  in  this  case 
it  was  found  that  d„  — »  0  as  u  — *  oo).  We  now  therefore  consider  the  homogeneous  system 
for  (3.34),  (3.35  a,b)  in  the  limit  as  v  — >  oc.  This  leads  to  a  system  of  difference  equations, 
and  we  therefore  surmise  that  the  coefficients  take  on  the  following  form  as  v  — ►  oc 

c„  ~  CqK  ",  (4.1) 

dv  «  D0K\  (4.2) 

whilst  because  of  the  nature  of  the  occurrence  of  the  we  must  have 

a„«  AQKve-ve\  (4.3) 

where  K  is  an  eigenvalue  whose  precise  value  is  extremely  important  since  it  is  this  value 
that  would  be  expected  to  determine  the  ultimate  downstream  behaviour  of  the  disturbance 
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field.  If  |  A' |  >  1,  then  using  the  generating  function  for  Bessel  functions  yields  the  following 


result 


~J2CoK'l'Jl'{z)  ~  C0e^K-l/K)\ 


This  then  predicts  exponential  growth  as  |rj  — >  oo  if  $l{z(K  —  1/A')}  >  0.  Therefore,  we 
proceed  to  the  actual  determination  of  K. 

Substituting  (4.1)-(4.3)  into  (3.34),  (3.35a,b),  and  discarding  inhomogeneous  terms, 
leads  to  the  following  three  equations  to  determine  K : 


cos(/3  -  6)  {[-eg  +  \(K  - 


—  - [2i  —  -  t~e°K  +  ee°K-l}AQ\ 
4  aa  J 


+  sin (/?  -  9)  |  ^P[e~0oK  +  e0* K'X)A0 

+ip2~C«  -  (ur  -  u2)[-—  +  Uh'  -  ^A'-’pol  =  0, 

U2  (is  2  2  J 


-  -[2*—  -  Ke~6°  +  K~1ee°]Ao 
4  a, 


-  cos(p  -  9)  (  l-P\Ke-do  +  A'-V°]A)  +  iP2~Co)\ 

14  U2  j  J  (4.6a,  b) 

+  {b22]  M0  ~  0)  {fi,  [~l^~  +  \i{K  -  A'"1)]  +  iP2as]Do} 

=  0. 

This  system  may  conveniently  be  expressed  in  the  form  of  a  sixth-order  linear  gener¬ 
alized  eigenvalue  problem  by  writing 

-do  =  A  i4p, 

Ch  =  A'Co,  (4.7) 

Do  =  A  Dq  , 
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and  solving  for  the  eigenvector  {.4o ,  .40,  Co,  Co,  ,  D0 } .  This  was  accomplished  numeri¬ 
cally  and  yielded  the  following  general  trends:  (i)  Of  the  six  roots  for  A  of  the  eigenvalue 
problem,  four  were  imaginary,  and  the  remaining  two  were  complex  twins  K  =  /v0  anti 
K  =  — c.c.{A’o}  (indeed,  the  nature  of  these  roots  can  be  confirmed  by  inspection  of  the 
system  (4.5),  (4.Ca,b);  (ii)  the  imaginary  roots  all  had  0  <  \K |  <  1:  (iii)  the  complex  pair  of 
roots  had  |A'|  >  1.  This  last  point  strongly  suggests  the  flow  is  structurally  unstable  since 
there  exists  homogeneous  solutions  for  o  and  F  which  grow  exponentially  as  the  distance 
from  the  wedge  tip  increases;  however,  further  detailed  analysis  given  later  shows  this  to 
be  untrue. 

Returning  now  to  consider  the  nature  of  the  imaginary  roots,  Figure  3  shows  the 
variation  of  these  roots  with  wedge  angle  6  at  M\  =  1.5,  3,  and  5.  In  addition  to  the 
trends  outlined  above,  the  roots  have  the  features  that  K  — >  0  as  6  — >  0.  whilst  as  0  — ►  6m, 

|  A' |  — *  1,  where  0m  is  the  limit  of  the  envelope  of  supersonic  flow  downstream  of  the  shock. 
This  then,  in  some  way  connects  with  the  results  of  previous  research  on  the  stability  of 
the  strong  shock  solution.  We  now  consider  the  nature  of  the  complex  roots.  Fortunately, 
it  is  possible  to  obtain  precise  analytic  information  on  these  roots  aided  by  the  observation 
that  the  eigenvectors  from  the  computation  reveal  that  D0  =  Dq  =  0  (to  within  machine 
precision).  If  we  use  this  result,  then  the  groupings  shown  in  (4.5),  (4.6)  would  suggest  that 
an  eigensolution  is  possible  if  the  following  three  equations  are  simultaneously  satisfied: 

+  l‘K  -  Uk~1]C0  -  \{2 1—  -  e-H'  4-  e8° K~']Ao  =  0,  (4.8) 

02  /  l  4  as 

l[e~doK  +  e^K-'jAo  +  ip^Co  =  0,  (4.9) 

4  U  2 

2*  -  —e~doK  +  —ee°K~'  =  0.  (4.10) 

a9  as 

We  find  this  is  indeed  possible  and  the  solution  is 

A'  =  '—{iaa±[Ul-ai}]'2}.  (4.11) 

Co  =  ±iAo/2.  (4.12) 
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Tilt'  value  of  K ,  given  by  (4.11)  has  | AT |  =  cs°  >  1,  suggesting  that  the  flow  is  always  struc¬ 
turally  unstable  in  the  sense  that  disturbances  will  always  grow  downstream.  However, 
it  turns  out  that  in  fact  this  homogeneous  solution  represents  a  structurally  stable  flow. 
This  stability  is  demonstrated  by  defining  the  following  functions  which  arise  naturally  as 
a  result  of  the  analysis  above 

E  =  C0eUJ~~~T^U^)  Y,  (4.13) 

t/=  —  oo 


together  with 


.  lu/t-oXo 


cf>  =  A0t  Y  {Kee'-6aY  J„{kr)-  (4.14) 

U—  — ‘OO 

E  and  <j>  are  both  harmonic  and  related  through  the  Cauchy-Riemann  equations,  with 


-  ~E; 


V2  ' 


<t>«  2  =  E 


*2- 


(4.15) 


As  a  direct  consequence  of  (4.15),  the  exponential  growth  of  the  two  components  to  this 
homogeneous  solution  cancel  each  other  when  computing  physical  quantities  such  as  the 
pressure  p  or  the  velocity  (u,5),  and  therefore,  the  flow  calculated  from  the  summations 
(3.26)-( 3.28)  is  bounded  for  all  values  of  (x2,y2). 

Although,  the  growing  homogeneous  solution  for  the  coefficients  does  not  lead  to  a 
physically  unbounded  solution,  due  to  exact  cancellations  when  calculating  physical  quan¬ 
tities,  their  presence  greatly  complicates  the  computation  of  the  summations  (3.26 )-( 3.28 ) 
for  calculating  the  flows  since  the  relevant  physical  quantities  result  from  the  remainder  of 
summing  three  or  more,  numerically  very  large  quantities,  to  yield  a  much  smaller  quantity. 
The  numerical  difficulties  in  obtaining  accurate  solutions  associated  with  this  approach  are 
obvious.  We  note  (again)  that  all  results  presented  have  been  carefully  checked,  to  account 
for  this  phenomenon. 

The  final  question  that  arises  is  the  nature  of  the  downstream  response  to  sustained 
acoustic  waves,  entropy  waves,  or  vorticity  waves.  We  may  expect  that  the  solution  down¬ 
stream  will  remain  oscillatory,  with  no  amplitude  decay,  although  this  is  not  always  the 
case.  The  different  downstream  behaviours  of  the  acoustic  disturbances  can  be  understood 
by  firstly  examining  the  solution  in  the  absence  of  the  wedge  (see  McKenzie  A'  Westphal 
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1968).  Since  the  frequency  of  the  disturbance  downstream  of  the  shock  is  fixed  by  the  fre¬ 
quency  of  the  disturbance  upstream  of  the  shock  and  the  wavenumber  of  the  disturbance 
tangential  to  the  shock  is  also  fixed  by  the  deflection  of  the  shock  produced  by  the  upstream 
disturbance,  there  is  a  finite  range  of  wavenumbers  tangential  to  the  shock  that  exclude  the 
existence  of  plane  acoustic  waves  which  propagate  downstream.  If  the  imposed  tangential 
wavenumber  lies  within  this  range,  then  the  acoustic  disturbance  behind  the  shock  (in  the 
absence  of  the  wedge)  exhibits  exponential  decay  in  a  direction  normal  to  the  shock.  If  the 
tangential  wavenumber  lies  outside  of  this  range,  then  a  plane  wave  propagates  at  an  angle 
0p  to  the  shock.  This  angle  and  the  range  of  excluded  tangential  wavenumbers  depend 
upon  the  relative  orientation  of  the  shock  and  the  upstream  base-flow  and  disturbance 
conditions. 

The  presence  of  the  wedge  changes  the  nature  of  these  solutions  drastically.  In  the 
case  of  the  tangential  wavenumber  being  within  the  excluded  range,  an  acoustic  field  is 
generated  that  decays  algebraically  (rather  than  exponentially)  owing  to  the  requirement 
that  the  shock  remain  attached  to  the  apex  and  the  requirement  that  disturbances  have 
zero  normal  velocity  at  the  wedge  surface.  If  the  tangential  wavenumber  is  outside  of 
the  excluded  range,  then  the  solution  in  the  presence  of  the  wedge  exhibits  two  types  of 
behaviour,  dependent  upon  the  angle  of  propagation  of  an  acoustic  disturbance  in  the 
absence  of  the  wedge.  If  the  angle  of  propagation  0p  is  greater  than  the  angle  between  the 
shock  and  the  wedge  0  —  0 ,  then  the  pressure  disturbances  generated  at  the  shock  would 
intersect  with  and  reflect  from  the  wedge  surface.  Thus,  there  would  be  a  non-decaying 
pressure  field  everywhere  between  the  shock  and  the  wedge  surface  with  a  superimposed 
algebraically  decaying  pressure  field  owing  to  the  requirement  of  an  attached  shock  and 
zero  velocity  normal  to  the  wedge  surface.  If  the  angle  of  propagation  0p  is  less  than  the 
angle  between  the  shock  and  the  wedge  surface,  then  the  pressure  field  is  divided  into  two 
regions  by  a  ray  emanating  from  the  apex  and  parallel  to  the  direction  of  propagation 
of  the  pressure  disturbance  in  the  absence  of  the  wedge.  Between  the  shock  and  this  ray 
there  is  a  non-decaying  pressure  field,  with  the  aforementioned  superimposed  non-decaying 
component.  Between  the  ray  and  the  wedge  surface  the  pressure  field  is  algebraically 
decaying  and  there  is  a  component  with  sustained  oscillations.  This  is  in  contrast  to  the 
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solution  in  the  absence  of  the  wedge  and  is  explained  by  the  presence  of  the  wedge  that 
serves  to  terminate  the  shock  (which  is  assumed  of  infinite  extent  in  both  directions  in  the 


absence  of  the  wedge),  and  there  are  no  pressure  disturbances  generated  by  the  shock  to 
fill  this  region.  We  will  refer  to  this  region  as  the  shadow  region. 

The  above  mentioned  downstream  limits  can  be  mathematically  understood  by  again 
considering  the  behaviour  of  the  coefficients  and  dv  as  u  — +  oc  when  the  inho¬ 
mogeneous  terms  i?L2\and  are  included.  For  this,  the  nature  of  9  in  (3.44)  is 

important.  If  9  is  real,  then  the  wavenumber  of  the  imposed  disturbances  lies  within  the 
excluded  range  derived  from  the  theory  of  McKenzie  and  Westphal  (1968),  the  coefficients 
within  the  summation  (3.43)  neither  grow  nor  decay,  and  the  combination  of  coefficients 
within  the  bracket  in  equation  (3.26)  all  decay  at  a  rate  proportional  to  e~l,Ba .  Thus  when 
reconstructing  the  pressure  along  any  ray  92  —  constant ,  algebraic  decay  is  found  since 
92  must  lie  on  the  interval  (0,#o)  with  zero  corresponding  to  the  wedge  surface  and  90 
corresponding  to  the  shock.  For  92  =  90  the  bracketed  combination  of  coefficients  does 
not  decay  leading  to  the  appropriate  oscillatory  behaviour  necessary  to  match  conditions 
at  the  shock.  If  9  is  complex,  then  the  condition 


(ori  cot  ft  +  o2)(cos  9  +  Asm#)  + 


iq-a] 


(4.17) 


is  satisfied.  We  also  note  that  for  equation  (3.43)  to  be  satisfied  then  sin#  must  remain 
real,  which  requires  that 


Notice  #,  is  related  to  the  angle  between  the  direction  of  the  shock  and  the  direction  of  the 
acoustic  disturbance  9p  through  tanh  9X  —  ft  tan  9p .  For  imaginary  values  of  9  the  coeffi¬ 
cients  of  (3.43)  behave  as  (±ietf’)*',  as  v  — >  oo,  which  increase  exponentially  in  magnitude 
as  do  the  corresponding  coefficients  R\}  \  R^\  and  R^\  If  9,  is  greater  than  #0,  then  the 
generated  or  transmitted  acoustic  waves  intersect  the  wedge  surface.  The  physical  conse¬ 
quence  of  this  intersection  is  seen  by  considering  the  bracketed  combination  of  coefficients 
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in  (3.26)  along  with  the  factor  cosh( vBz)-  The  product  grows  exponentially  and  alternates 
between  real  and  imaginary.  The  relation  (4.4)  shows  that  there  is  a  sustained  oscillatory 
solution  of  the  pressure  for  all  x-i  and  y2-  For  0  <  0,  <  the  aforementioned  product 
of  coefficients  is  exponentially  growing  when  02  >  0 1  which  shows  that  there  is  a  region  of 
sustained  oscillatory  behaviour  of  the  pressure  field  near  the  shock,  and  the  aforementioned 
product  of  coefficients  is  exponentially  decaying  for  O2  <  0t  which  shows  that  there  is  a 
region  of  algebraic  decay  of  the  pressure  field  near  the  wedge  surface  and  that  there  are 
no  sustained  oscillations  near  the  surface.  This  overall  behaviour  of  the  pressure  field  has 
been  verified  by  direct  computation  of  the  summation  (3.26)  along  rays  62  =  constant  for 
various  cases. 

With  the  previous  analysis  in  mind  and  restating  that  one  goal  of  this  analysis  is  to 
discover  how  these  solutions  might  affect  the  boundary  layer  receptivity  problem,  we  focus 
attention  on  the  pressure  evaluated  along  the  wedge  which  generally  may  naturally  be 
decomposed  into  four  parts 


e  lu,lp  =  P0e  '2  a» 


+  P\x2  1/2e  ( 
+  P2x~1/2e~{ 


iu/  x  2 
U 2- a 


->t) 


iu/x  2 


-*f) 


4-  Q{x 2). 


(4.20) 


The  first  term  on  the  right-hand-side  is  non-zero  only  when  0,  >  0o  and  is  evaluated 
numerically  by  considering  the  large  v  behaviour  the  coefficients  {«„ }.  This  term  is  con¬ 
verted  to  an  infinite  sum  of  Bessel  functions  and  then  subtracted  from  the  left-hand-side 
where  the  pressure  has  the  representation  given  by  the  summation  (3.26).  The  Fourier 
Sine  Transform  of  the  resulting  difference  is  evaluated  and  the  strengths  of  the  inverse 
square  root  singularities  provide  the  amplitudes  P\  and  P2  for  the  second  and  third  terms. 
Physically,  the  first  term  represents  the  transmitted  acoustic  wave  (or  generated  acoustic 
wave  if  the  imposed  disturbance  is  an  entropy  or  a  vorticity  wave)  when  it  intersects  the 
wedge  surface,  and  the  second  and  third  term  are  the  far-downstream  behaviour  of  slow 
and  fast  moving  acoustic  waves  associated  with  reflections  off  the  wedge  surface  and  regen¬ 
eration  of  new  acoustic  waves  at  the  shock,  as  well  as  new  acoustic  waves  generated  by  the 
requirement  that  the  shock  remains  attached  to  the  wedge  tip.  The  first  three  terms  are 
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tilt*  leading  terms  in  the  limit  x2  — >  oc,  and  the  fourth  term  decays  faster  than  J'7*  ‘  and 
is  considered  to  he  the  local  field.  The  Fourier  spectrum  of  the  fourth  term  is  important 
since  it  represents  the  generation  of  all  wavelengths  by  the  response  of  the  attached  shock 
to  a  single  wavelength  input  ;  the  spectrum  of  the  local  field  is  hounded  for  all  wavenumbers 
since  the  local  field  decays  relatively  quickly  as  x2  — >  oc. 

\{0,  <  then  the  far-downstream  behaviour  of  the  pressure  along  the  surface  changes. 
The  wavenumber  of  the  imposed  disturbance  is  outside  of  the  allowable  wavenumbers  for 
propagating  plane  waves  or  the  generated  or  transmitted  plane- wave  pressure  disturbance 
does  not  intersect  the  wedge,  and  it  is  found  that  P0  is  zero.  The  amplitudes  Pj  and  P2 
are  11011-zero  and  therefore  the  pressure  has  algebraic  decay  far  downstream. 

The  vorticity,  found  by  computing  =  V2P,  can  be  similarly  decomposed  into  four 
parts  if  6,  >  0 

_  A  y  2 

_  e~‘'2<l'f-*?){ft0(d*cosh  0,1,2 
+  Q2i/2-,/2e*^+‘T}+9(j/2), 

where  we  have  assumed  that  the  calculations  take  place  far  enough  downstream  in  x2  such 
that  the  behaviour  represented  by  (4.21)  is  an  accurate  expansion  valid  for  large  y2 .  The 
first  term  represents  a  generated  or  transmitted  vorticity  wave  and  the  next  two  terms 
represent  additional  vorticity  generated  by  the  moving  shock  as  a  result  of  requiring  that 
the  shock  remain  attached  to  the  wedge  apex  and  zero  normal  velocity  on  the  surface.  That 
is,  the  algebraically  decaying  component  of  the  pressure  field  generated  by  the  presence  of 
the  wedge  interacts  with  the  shock  wave  through  the  Rankine-Hugoniot  relations  which 
indicate  that  the  shock  must  be  distorted  by  this  component  <  f  the  pressure  field,  and 
hence  algebraically  decaying  vorticity  modes  represented  by  ill]  and  fi2  rre  generated 
by  the  curvature  of  the  shock.  A  numerical  subtlety  exists  when  trying  to  calculate  the 
decomposition  (4.21)  when  0,  >  0.  The  coefficient  of  the  vorticity  depends  on  cancellations 
involving  five  very  large  terms  (for  the  pressure  decomposition  there  were  only  three  very 
large  terms);  thus,  it  is  not  possible  to  accurately  determine  flo  for  some  cases.  If  0,  =  0. 
then  it  appears  to  be  not  possible  to  calculate  a  decomposition  of  the  vorticity  into  the  form 
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given  in  (4.21).  The  coefficients  of  the  vorticity  expansion  neither  grow  nor  decay  making 
the  numerical  determination  of  difficult  if  not  impossible.  In  the  absence  of  the  wedge, 
a  single  wavelength  vorticity  mode  would  still  exist  for  this  case,  and  the  pressure  field 
would  be  exponentially  decaying.  With  the  wedge  present,  there  exists  a  small  algebraic 
component  of  the  pressure  field. 

The  results  relating  to  the  decompositions  given  in  (4.20)  and  (4.21 )  are  given  in  the 
next  section. 


5.  Results 

The  response  of  the  flow  downstream  of  the  shock  is  generally  complicated,  although 
the  analysis  of  the  previous  section  does  enable  us  to  quantify  the  response,  at  least  in 
the  far  downstream  limit.  To  this  end,  we  find  it  useful  to  consider  the  sine- transform  of 
Q(x 2 ),  that  is 

Q*  =  I  (J(x2)sin(p:r2)<i:r2.  (5.1) 

Jo 

Fortunately,  this  quadrature  can  be  carried  out  semi-analytically  using  the  following  sine- 
transform 


i; 


sm(/ix2)e,?I?  Jv(kx2)dx,2  = 


i 


(p  +  q )2  -T  i(n  +  q) 


\ 


2  ikl 


/F-(^  +  9)2 

j  |  \A2  -  &  -  q y  -  *(p  -  ?)  j 


(5.2) 


2  ikv 


\A2  -  0*  -  ?>2 


The  inverse  square-root  singularities  lead  to  the  second  and  third  terms  of  (4.20)  and 
the  evaluation  of  (5.1)  involves  a  summation  of  terms  derived  from  (5.2)  with  coefficients 
determined  numerically.  The  results  of  this  process  are  shown  graphically  in  Figures  4-8 
where  3£{Q*}  is  denoted  by  solid  lines,  whilst  Cr{Q*}  is  denoted  by  broken  lines. 

The  spectrum  Q*  for  the  input  conditions  M\  =  5,  ot\  —  1,  o2  =  0  and  6  =  25° 
(slow  mode)  is  shown  in  Figure  4.  For  this  set  of  input  conditions,  inequality  (4.17)  is  not 
satisfied  and  therefore  Pq  in  (4.20)  is  zero.  This  spectrum  is  bounded  for  all  wavenumbers 
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Imr  has  square-root  behaviour  at  //  =  .  the  wavenumber  location-.  of  the  slow  ami 

fast  acoustic  disturi  mces  behind  the  shock.  This  square  root  behaviour  remain*-  alter 
subtracting  the  inverse  square-root  singularity  to  form  the  second  and  third  terms  of  ( 4.20 ). 
For  wavenumbers  between  these  two  pemts.  a  plan**  wave  with  the  fixed  rim*'  frequency  w 
would  have  to  travel  faster  than  the  sound  speed  relative  to  the  mean  How.  and  so  th.  ‘Sr 
represent  pressure  disturbances  that  decay  exponentially  in  x  >-  The  local  Held  is  comprised 
of  a  Fourier  integral  over  both  classes  of  disturbance.  It  should  be  noted  that  the  large 
wavenumber  disturbances  are  “slow  modes"  whilst  the  small  wavenumber  disturbances  air 
"fast  modes”.  The  importance  of  Q*  is  that  it  demonstrates  that  far  downstream  the  energy 
density  of  the  acoustic  :  ponse  is  spread  over  all  wavenumbers  and  is  not  concentrated 
at  an  isolated  wavenumber  as  is  the  energy  density  of  tin*  acoustic  mode  in  t lit'  absence  of 
the  wedge.  In  particular,  the  energy  density  is  significant  at  large  wavenumbers,  i.r.  small 
wavelengths,  and  these  disturbances  generated  by  tin*  shock  body  intera  .  *on  might  lead 
to  receptivity  as  they  interact  with  a  growing  boundary  layer  on  the  wedge  surfacm 

Figure  5  shows  the  spectrum  for  the  fast  mode  with  tin*  same  input  conditions  as  in 
Figured.  For  the  fast  mode,  inequality  (4.17)  is  satisfied  so  that  Pu  is  non-zero.  To  indicate 
tin*  presence  of  this  transmitted  plant*  wave,  a  circle  has  been  added  to  Figure  5  at  the 
wavenumber  location  of  the  transmitted  wave.  It  should  be  noted  that  this  wave  number 
location  is  not  quite  at  the  same  wavenumber  as  the  singularity  of  the  Fourier  transform 
which  leads  to  the  second  and  third  terms  of  the  expansion.  Figure  4  and  Figure  5  appear 
to  be  quantitatively  similar,  although  a  shift  to  smaller  wavenumbers  is  forced  by  the 
changed  frequency  of  the  imposed  disturbance,  a  result  of  choosing  the  opposite  sign  in 
equation  (2.8). 

In  Figure  G.  the  spectrum  of  a  fast  mode  with  larger  incoming  Mach  number  is  shown. 
The  input  conditions  are  M\  =  10.  tq  =1,  o-2  =  0  and  0  =  25°  (slow  mode).  This 
set  of  input  conditions  also  falls  within  the  envelope  of  flows  which  lead  to  sustained 
downstream  acoustic  waves.  It  is  seen  that  the  energy  density  of  the  local  field  has  a 
strong  Mach  number  dependence,  as  do  the  amplitudes  of  the  first  three  terms  of  (4.20) 
which  will  be  explored  in  Figures  9-12.  Figures  7  and  8  show  the  spectra  for  the  local 
pressure  field  Q* ,  in  response  to  an  incoming  entropy  disturbance  (2.13-2.1G)  and  to  an 
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incoming  vorticity  disturbance  (2.9-2.12).  The  input  conditions  arc  the  same  as  in  Figur> 

4  (except  for  Figure  8  for  which  we  have  taken  (ll  —  ()  and  o _>  =  1  ).  There  is  a  qualitative 
difference  between  the  local  field  response  to  acoustic  disturbances  and  tin  local  field 
response  to  these  non-acoustic  inputs.  The  spectra  in  Figures  7  and  S  are  considerably  more 
peaked  near  the  fast-mode  singularity  with  the  response  near  the  slow-mode  singularity  not 
nearly  as  strong.  For  the  re:  poase  to  both  the  fast  and  slow  incoming  acoustic  modes,  the 
spectra  nem  both  singularities  were  broad  with  the  strength  of  the  slow-mode  singularity  as 
strong  or  stronger  than  the  fast  mode  singularity.  Perhaps  this  is  a  result  of  the  incoming 
entropy  disturbance  having  zero  pressure  and  velocity  components,  the  incoming  vorticity 
wave  having  zero  pressure  and  density  components,  whilst  all  components  of  the  incoming 
acoustic  disturbances  are  non-zero,  thus  leading  to  a  much  more  complicated  interaction 
with  the  attached  shock  (especially  near  tin*  tip  of  the  wedge). 

The  amplitude  of  the  transmitted  or  generated  (sustained)  acoustic  disturbance  TV 
the  amplitude  of  the  algebraically  decaying  slow  modi'  P\ .  and  the  amplitude  of  the  alge¬ 
braically  decaying  fast  mode  P2  are  shown  graphically  in  Figures  9-12.  It  can  be  shown 
that  To  and  P2  are  real  quantities  whilst  P\  is  strictly  imaginary.  9?{  Pi }  is  denoted  as  a 
solid  line,  is  denoted  as  a  line  with  small  dashes  and  4?{p>}  >s  denoted  as  a  line 

with  large  dashes. 

Figure  9a  shows  the  response  for  the  an  impinging  slow  acoustic  mode.  Over  a  large 
range  of  Mach  numbers  of  the  incoming  flow  M\ ,  there  exists  no  downstream  sustained 
acoustic  waves,  i.e  Po  —  0;  however,  both  of  the  algebraically  decaying  modes  are  non-zero, 
and  the  amplitude  s  increase  considerably  as  the  Mach  number  increases.  As  mentioned  in 
tin'  introduction,  this  is  significant  since  in  the  absence  of  the  wedge  the  theory  (as  verified 
by  numerical  computations)  uredicts  that  the  acoustic  response  will  be  exponentially  de¬ 
caying,  but  by  requiring  that  the  shock  remains  attached  to  the  ’edge  tip.  it  is  seen  that 
the  response  decays  like  the  inverse  square  roof  (together  with  a  local  field  which  also  does 
not  decay  exponentially).  Figure  9b  shows  the  details  of  the  transition  region  (in  Mach 
number)  between  the  decaying  pressure  response  and  the  sustained  pressure  response.  The 
amplitudes  of  P\  and  P>  jump  dramatical’  when  P()  is  non-zero.  Physically,  this  implies 
that  then'  are  two  parts  to  the  algebraically  decaying  modes.  Tilt'  first  part  is  a  result  of 
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tlu'  requirement  that  the  shock  remain  attached  as  seen  when  Po  =  0;  the  second  part  is 
a  direct  result  of  the  sustained  acoustic  disturbance  reflecting  off  the  wedge  surface  and 
interacting  with  the  shock  to  generate  more  acoustic  disturbances.  As  A/]  decreases,  the 
critical  wedge  angle  (above  which  there  is  no  mean  flow  solution)  decreases  until  it  reaches 
the  fixed  wedge  angle  of  Figures  9-12  ( i.e. .  9  =  25°).  The  amplitude  of  the  transmitted 
wave  remains  finite  in  this  limit;  however  P\  and  Pi  increase  sharply  in  magnitude,  this  im¬ 
plies  that  the  solution  is  becoming  increasingly  unstable  as  the  critical  angle  is  approached 
(as  would  be  expected).  In  Figures  10-12,  the  response  to  an  incoming  fast  acoustic  wave, 
to  an  incoming  entropy  wave  and  to  an  incoming  vorticity  wave,  respectively,  are  shown; 
in  all  cases  9  =  25°,  M\  =  5,  aj  =1.  and  cv2  =  0  except  in  the  case  of  the  vorticity  mode 
(Figure  12)  where  0-2  =  1  and  oq  =0.  There  is  a  sustained  acoustic  disturbance  for  these 
three  cases,  and  correspondingly,  the  amplitudes  of  P\  and  P2  are  relatively  large  since  as 
mentioned  in  the  discussion  of  Figure  9  these  two  modes  have  components  related  to  the 
reflection  of  the  acoustic  waves  from  the  wedge  surface.  All  three  figures  demonstrate  the 
same  behaviour  as  Figure  9  when  the  critical  angle  is  approached  -  the  amplitude  of  the 
sustained  mode  remains  finite  whilst  the  other  modes  grow.  For  the  responses  to  the  im¬ 
pinging  acoustic  disturbance  (Figure  10)  and  to  the  impinging  entropy  disturbance  (Figure 
11),  the  algebraically  decaying  fast-mode  is  considerably  larger  than  the  slow-mode  with 
the  relative  importance  of  the  slow-mode  being  greater  for  the  response  to  the  entropy 
wave.  Comparing  the  spectra  shown  in  Figures  5  and  Figure  7,  the  energy  of  the  :Aow 
mode  is  seen  to  be  in  the  local  fields  with  relatively  more  energy  for  the  slow-mode  in  the 
response  to  the  impinging  acoustic  wave.  I11  the  response  to  the  incoming  vorticity  distur¬ 
bance  (Figure  12),  the  amplitude  of  the  slow  mode  P\  is  much  larger  than  the  amplitude 
of  the  fast-mode  P2.  Comparing  with  Figure  8,  this  is  consistent  with  the  relative  absence 
of  energy  in  the  slow-mode  of  the  local  field. 

For  completeness.  Figure  13  and  Figure  14  show  the  result  of  a  decomposition  of  the 
vorticity  into  the  form  given  by  (4.21).  Figure  13  represents  a  response  to  an  incoming 
fast  acoustic  wave  with  the  same  conditions  as  Figure  10,  but  with  a  wedge  angle  of  only 
20°.  The  change  of  wedge  angle  is  due  to  the  numerical  difficulties  of  calculating  al¬ 
ready  mentioned.  Figure  14  shows  the  response  to  an  incoming  entropy  wave.  General 
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conclusions  about  the  decomposition  are  difficult  to  obtain  (again  related  to  the  numerical 
difficulties);  however,  we  do  note  the  relative  insensitivity  to  the  Mach  number  of  the  eoef 
ficients  except  near  the  critical  angle  where  one  of  the  algebraically  decaying  components 
becomes  quite  large.  Again,  this  probably  relates  directly  to  the  stability  of  the  entire 
solution  as  the  critical  angle  is  approached. 


6.  Conclusions 

We  have  considered  the  interaction  between  the  shock  wave  attached  to  the  wedge  and 
freest  ream  disturbances.  In  order  to  quantify  the  results,  the  freestream  disturbances  were 
taken  to  be  fixed-frequency,  single-wavelength,  plane-wave  disturbances  of  acoustic  type, 
the  vorticity  type,  or  entropy  type.  The  behaviour  of  the  flow  field  downstream  of  the  shock 
in  the  shock- wedge  combination  was  compared  with  the  behaviour  of  the  downstream  flow 
in  the  absence  of  the  wedge.  After  rationalizing  a  number  of  mathematical  and  computa¬ 
tional  subtleties,  we  were  able  to  show  precisely  how  the  presence  of  the  wedge  affects  the 
interaction  between  a  shock  and  a  free  stream  disturbance.  We  specifically  focused  on  the 
pressure  disturbance  along  the  wedge  surface  as  the  analysis  of  these  disturbances  directly 
relates  to  boundary-layer  receptivity  in  the  shock-wedge  configuration  and  we  also  studied 
the  vorticity  generated  by  the  presence  of  the  wedge  as  this  analysis  relates  to  enhanced 
mixing  by  the  vorticity  production  of  a  moving  shock. 

The  most  significant  differences  in  the  interaction  produced  by  the  presence  of  the 
wedge  are  the  presence  of  an  acoustic  and  vorticity  field  which  decays  algebraically  in 
addition  to  sustained  oscillations  produced  by  the  incoming  disturbance.  It  has  been 
shown  that  this  field  is  comprised  of  two  parts:  one  produced  by  the  requirement  that  the 
shock  remain  attached  to  the  apex  of  the  wedge  and  the  other  produced  by  reflections  of 
acoustic  waves  from  the  wedge  surface. 
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p  for  case  1  (acoustic  mode) 
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M 

Fig.9a  Po,  Pi,  P2  upstream  slow  acoustic  mode,  0  =  25°,  aj  =  1, 
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p  for  case  2  (acoustic  mode) 


Fig.  10  P0,  Pi,  P2  for  upstream  fast  acoustic  mode,  8  =  25°,  c*i  =  1,  a2  =  0. 
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Fig. 12  Pot  Pi,  P2  for  upstream  vorticity  mode,  9  =  25°,  ai  —  1,  <22  —  0- 
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vorticity  coefficients  for  case  2  (acoustic  wave), 13  =  20 
3  r 
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vorticity  coefficient 


vorticity  coefficients  for  cose  3  (entropy  wove) 


1 

Fig. 14  fio,  ,  ^2  for  upstream  entropy  mode,  6  —  25°,  <*i  =  1,  ct2  —  0 
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